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Abstract 



We present a framework for discussing LES equations with nonlinear dispersion. In this 
framework, we discuss the properties of the nonlinearly dispersive Navier-Stokes-alpha 
(NS— a) model of incompressible fluid turbulence — also called the viscous Camassa- 
Holm equations in the literature — in comparison with the corresponding properties of 
large eddy simulation (LES) equations obtained via the approximate-inverse approach. 
In this comparison, we identify the spatially filtered NS— a equations with a class of 
generalized LES similarity models. Applying a certain approximate inverse to this class 
of LES models restores the Kelvin circulation theorem for the defiltered velocity and 
shows that the NS— a model describes the dynamics of the defiltered velocity for this 
class of generalized LES similarity models. We also show that the subgrid scale forces 
in the NS— a model transform covariantly under Galilean transformations and under 
a change to a uniformly rotating reference frame. Finally, we discuss in the spectral 
formulation how the NS— a model retains the local interactions among the large scales ; 
retains the nonlocal sweeping effects of large scales on small scales ; yet attenuates the 
local interactions of the small scales amongst themselves. 

1 Introduction 

The effects of subgrid-scale (SGS) fluid motions occurring below the available resolution of 
numerical simulations must be modeled. One way of modeling these effects is simply to 
discard the energy that reaches such subgrid scales. This is clearly unacceptable, though, 
and many creative alternatives have been offered. A prominent example is the large eddy 
simulation (LES) approach, see, e.g., IQ - 0- The LES approach is based on applying a 
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spatial filter to the Navier-Stokes equations. The reduction of flow complexity and informa- 
tion content achieved in the LES approach depends on the characteristics of the filter that 
one uses, its type and width. In particular, the LES approach introduces a length scale into 
the description of fluid dynamics, namely, the width of the filter used. The LES approach 
is conceptually different from the Reynolds Averaged Navier-Stokes (or, RANS) approach, 
which is based on statistical arguments and exact ensemble averages, rather than spatial 
and temporal filtering. After filtering, however, just as in the RANS approach, one faces 
the classic turbulence closure problem: How to model the effects of the filtered-out subgrid 
scales in terms of the remaining resolved fields? In practice, this problem is compounded 
by the requirement that the equations be solved numerically, thereby introducing further 
approximations. 

A turbulence modeling scheme — called here the Navier-Stokes-alpha model, or NS— a 
model (also called the viscous Camassa-Holm equations in - [^) — was recently in- 
troduced. This model imposes an energy "penalty" inhibiting the creation of smaller and 
smaller excitations below a certain length scale (denoted alpha), but still allows for nonlinear 
sweeping of the smaller scales by the larger ones. This energy penalty implies a nonlin- 
early dispersive modification of the Navier-Stokes equations. The alpha-modification 
appears in the nonlinear convection term, it depends on length scale and we emphasize that 
it is dispersive, not dissipative. This modification causes the translational kinetic energy 
wavenumber spectrum of the NS— a model to roll off rapidly below the length scale alpha as 
in three dimensions, instead of continuing to follow the slower Kolmogorov scaling law, 
k~^^^, thereby shortening its inertial range and making it more computable |]10 . 



This roll-off in the energy spectrum means that the inertial range is curtailed and the 
length scale ia at which viscous dissipation takes over in the NS— a model is larger than for 
the Navier-Stokes equations with the same viscosity. Hence, for a given driving force and 
viscosity, the number of active degrees of freedom N^^j for the NS— a model is sma//er than 
for Navier-Stokes. A rigorous analytical estimate was derived in P] for the fractal dimension 
Dfrac of the global attractor for the NS— a model. Namely, 

Dfrac<{NZfr/' and N^f ^ {L/e^f c:^ ^Re'/\ (1) 

where L is the integral scale (or domain size), ia is the end of the NS— a inertial range 
and Re = L^^^el/^/u is the Reynolds number (with NS— a energy dissipation rate Sa and 
viscosity u). The Kolmogorov dissipation length scale is ixo = (i^^/^aY^^- The three length 
scales a > ia > ixo are related by 



«4o or -c^j^. (2) 



Thus, these length scales a > ia > ixo stand in relation to each other in the same way as 
the length scales {K^/e'^Y^^ > {Ku/eY^'^ > {u^/eY^'^ in a K-epsilon model of Navier-Stokes 
turbulence. 
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The number of degrees of freedom for a corresponding Navier-Stokes flow with the same 
viscosity and energy dissipation rate is 

NZ'f ^ {L/iKof - Re'" , (3) 

The imphcation of these estimates of degrees of freedom for direct numerical simulations that 
access a significant number of them using the NS— a model is an increase in computational 
speed relative to Navier-Stokes of 



(f)'"=^(i)'-- 



Thus, if a tends to a constant value, say L/lOO, when the Reynolds number increases - as 
found in p - by comparing steady NS— a solutions with experimental data for turbulent 
flows in pipes and channels - then one could expect to obtain a substantial increase in 
computability by using the NS— a model at high Reynolds numbers. An early indication of 
the reliability of using these estimates to gain a relative increase in computational speed in 
direct numerical simulations of homogeneous turbulence in a periodic domain is given in 0. 

In this paper, we introduce the NS— a model as one of a class of equations obtained 
by filtering the Kelvin circulation theorem for Navier-Stokes. These equations possess the 



vortex stretching and transport properties emphasized in |T^. We then compare a spatially 
filtered version of the NS— a equations with a generalization of the LES similarity models. 
We also calculate the kinematic tranformation properties of the NS— a model under changes 
of reference frames and discuss its nonlinear spectral-mode dynamics. 



2 Kelvin-filtered turbulence models 

Following [|TU|, the Navier-Stokes-alpha (NS— a) model may be introduced by starting from 
Kelvin's circulation theorem for the original Navier-Stokes (NS) equations, 

— + V ■ Vv + Vp = t^V V + f , with V ■ V = , (5) 
at 

for an appropriate forcing f and constant kinematic viscosity u. These equations satisfy 
Kelvin's circulation theorem, 

4 / v-dx= <f (z/V^v + f ) ■ c/x , (6) 

at J-yiv) J-y{v) 

for a fluid loop 7(v) that moves with velocity v(x, t), the Eulerian fluid velocity. 



Kelvin-filtering the Navier-Stokes equations. The equations for the NS— a model 
may be introduced by modifying the NS Kelvin circulation theorem to integrate around a 
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loop 7(v) that moves with a spatially filtered Eulerian fluid velocity given by v = g*v, 
where * denotes the convolution, 

ir = g*^r = | (^(x - x' ) v(x' ) d^^' = L,v(x) . (7) 

The "inverse" operation is denoted 

v = Ov = L^-iv, (8) 

thereby defining an operator O whose Green's function is the filter g and which we shall 
assume is positive, symmetric, isotropic, translation-invariant and time-independent. Under 
these assumptions the quantity (kinetic energy) 

E = ^ J V ■ V d^x = ^ J V • 5( * V d^x ^ \ j ^ ' ' (9) 

defines a norm. 

We modify the Navier-Stokes equations (|^) by replacing in their Kelvin's circulation 
theorem (^) the loop 7(v) with another loop 7(v) moving with the spatially filtered velocity, 
V. Then we have, 

4- / V ■ rfx = / (z/VV + f ) ■ rfx . (10) 

dt J-y(ir) J7(v) 

After taking the time derivative inside the Kelvin loop integral moving with filtered velocity 
V and reconstructing the gradient of pressure, we find the Kelvin- filtered Navier-Stokes 
equation, 

— + v- Vv + Vv'^-v + Vp = z/V^v + f , (11) 

with 

V ■ V = , and v = Ov . (12) 

The transport velocity v(x, t) is the spatially filtered Eulerian fiuid velocity in Eq. (^. Note 
that the continuity equation is now imposed as V ■ v = V ■ ((? * v) = 0. If filtering commutes 
with differentiation, then (7 * (V ■ v) also vanishes and the transported velocity v is also 
divergenceless. The energy balance relation derived from the NS— a equation (p!T| ) is 

^ y" i V ■ V rf^s = J V ■ i d^x - u J tr(Vv^ • Vv) d^x , (13) 

where we have used incompressibility and have dropped any boundary terms that appear 
upon integrating by parts. 

Vortex transport and stretching. Let q = V x v be the vorticity of the unfiltered 
Eulerian fluid velocity. The curl of the Kelvin-filtered Navier-Stokes equation ([TT|) gives the 
vortex transport and stretching equation, 

^ + v-Vq-q- Vv = z/V^q + Vx f. (14) 
ot 
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We note that the coefficient Vv in the vortex stretching term is the gradient of the spa- 
tially filtered Eulerian velocity v. Thus, Kelvin-filtering tempers the vortex stretching in 
the modified Navier-Stokes equation (|ll]), while preserving the original form of the vortex dy- 
namics. This tempered vorticity stretching is reminicent of Leray's approach in regularizing 
the Navier-Stokes equations. Except for the term (Vv)""" ■ v, the Kelvin-ffitered Navier- 
Stokes equation ( pAf ) is otherwise quite similar to Leray's regularization of the Navier-Stokes 
equations proposed in 1934 |ll|. Extension of the Leray regularization to satisfy the Kelvin 



circulation theorem was cited as an outstanding problem in Gallavotti's review |T2|. Indeed, 
for the case that O is the Helmholtz operator, full advantage of this regularized vorticity 
stretching effect was taken in to prove the global existence and uniqueness of strong solu- 
tions for the three dimensional NS— a model. Paper also proves the estimate (|l|) for the 
fractal dimension of the global attractor for the NS— a model. 

The Navier-Stokes- alpha model emerges from the Kelvin-filtered Navier-Stokes equa- 
tion (IT^) when one chooses the operator O to be the Helmholtz operator, thereby introducing 
a constant a that has dimensions of length, 

= 1- a^V^ , with a = const . (15) 

In this case, the filtered and unffitered fluid velocities in Eq. (|D are related by 

v = (l-a2v2)v, (16) 

and the reconstructed pressure p is related to the hydrodynamic pressure P as 

p = P-l|vp-^|Vvr. (17) 

The original derivation of the ideal Euler-alpha model (the v = Q case of the NS— a 
model) was obtained by using the Euler-Poincare approach in |T^, [Q. The physical in- 



terpretations of V and v as the Eulerian and Lagrangian mean fluid velocities were given in 
The kinetic energy norm corresponding to (^ for the NS— a model is given by 



0. 



a 



2 



1 , . ,9 a , 

- vr H Vv 

2' ' 2 ' ' 



(fx. 



This kinetic energy is the sum of a translational kinetic energy based on the spatially ffitered 
Eulerian velocity v, and a gradient-velocity kinetic energy, multiplied by a^. By showing 
the global boundedness, in time, of the kinetic energy (0) one concludes that the coefficient 
Vv in the vortex stretching relation (0) is hounded in for the NS— a model. Thus, the 
second term in the kinetic energy (p!8D imposes an energy penalty for creating small scales. 
The spatial integral of | Vvp in the second term has the same dimensions as filtered enstrophy 
(the spatial integral of |V x vp, the squared filtered vorticity). For a domain with boundary, 
or when considering transformations to rotating frames, the spatial integral of || Vvp in the 
second term in the kinetic energy norm (|18|) should be replaced by the integral of trace (e-e). 
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where e is the strain rate tensor, e 
coordinates 



16 



(1/2) (Vv + Vv""") of the filtered velocity in Euclidean 
For the case we treat here, in Euclidean coordinates and in the absence of 



boundaries, all three of these norms are equivalent for an incompressible flow. 



3 Comparison of the NS— a model with LES equations 

The large eddy simulation (LES) equations for an incompressible flow are conventionally 
written as 

d _ d d ^ _ d 

-^Ui + —UiUj = -—P + U———Ui- —Tij (19) 
Ot OXj OXi OXjOXj OXj 

with -7r^ = 0. 20) 

OXi 

The overbar denotes spatial filtering which, for a quantity /(x), is defined by the integral 

7(x) = I G(x - x' ) /(x' ) dV = Lg/(x) , (21) 

where G is a given kernel function, not necessarily the same as g in Eq. The quantities 
in Eq. ([I9|) , Ui, P, and u are the velocity components, pressure, and the kinematic viscosity, 
respectively, and the effects of subgrid-scale quantities on the resolved velocity are described 
by a subgrid-scale (SGS) stress tensor 

Tij = v~Uj -UiUj, (22) 

which must be modeled in terms of the resolved quantities to close Eq. ([T9|). Note that 
the unfiltered velocity Ui is assumed to be a solution of the incompressible Navier-Stokes 
equation. 

In contrast, the NS— a equation ([TT|) with definition (0) and (p^Tf ) is written in index 
notation as 

dt * ^ dxj dxi dxi \ 2 2 dxi dxi J dxjdxj 

where Vi is the filtered quantity related to the NS— a velocity Vi through the equation 

y.={j.- y. =Vi- C?\j'^Vi (24) 

OXjOXj 

and a is the filter length scale. Both fields Vi and are incompressible since the 'hat' filter 
commutes with differentiation for a constant alpha. 

The NS— a equation ( p3| ) can be rearranged into a form similar to the LES equation ( |T9|) 

^y. + —y.y. = - — + ^, ^ —iiv- -v-)v- - a^^^\ (25) 

dt * dxj * ^ dxi dxjdxj * dxj \ ^ ^ ^ dxj dxi ) 
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Two interpretations of the NS— a velocity. Both the NS— a velocity Vi and the filtered 
NS— a velocity Vi can be represented numerically using a lower resolution than required to 
represent the full NS velocity for the same problem. This is the property that they share 
with the filtered velocity in traditional LES. Nevertheless, it must be recognized from the 
derivation in the previous section that there is no LES filtering procedure that would produce 
the NS— a equations from the full NS equations. Nevertheless because of the similar behavior 
both NS— a velocities can be formally interpreted as the traditional LES velocity. If the 
NS— a velocity f « in equation (p5D is interpreted as the filtered Navier-Stokes velocity Ui in 
the LES equation ([l9|) then the subgrid-stress tensor associated with the alpha model is 



where the last equality is obtained using Eq. (p4|) . This is a straightforward identification 
of terms. 

Alternatively, one may propose to interpret the filtered NS— a velocity Vi as the filtered 
Navier-Stokes velocity Wj. To obtain the corresponding SGS stress in this case we first filter 
the NS— a equation (|25|) with the 'hat' filter. Since this filtering commutes with differentia- 
tion one finds 



9 a _ dP 92 ^ d (r^ _x ^(dvkdvk , /OVA 
^^i + ^t-^t-j = - — + z/ Vi-— \iviVj - VivA - a [-^-^ VVj\/^Vi] \ , (27) 

Ot OXj OXi OXjOXj OXj \^ ' \OXj OXi J J 

which involves doubly filtered quantities. The corresponding SGS stress tensor is 



rlf = {viVj - ViVj) - a" (|^|^ + ^V^i)^) (28) 
= {viVj-ViV^)-a^^^^-a'(vj\/Hi-v,\/''v^ (29) 

where the last equality is obtained using Eq. (p^). We shall next consider the transformation 
properties of the NS— a model. In section |^ we shall discuss how the SGS stress tensor r/," 
for the filtered NS— a model is related to LES similarity models for NS. 



4 Transformation properties 

Galilean invariance. Because of the first term in Eq. (^), or equivalently the expression 
ViV^-Oj in the second part of (p6|), the SGS tensor rf^ is neither symmetric nor Galilean 
invariant, in contrast to the ordinary SGS stress tensor ([2^). However, the SGS force which 
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is the divergence of the SGS stress tensor is Gahlean invariant for both representations of 
T^j. Therefore, the dynamics of the alpha model is Galilean invariant. Similarly, the SGS 
stress tensor r/^ is not symmetric and lacks the Galilean invariance because of the term 
Vi'V'^Vj but the Galilean invariance is restored in its divergence since Vj is incompressible. 

Rotating frame. It is of interest to consider a transformation of the SGS stress tensor 
to a frame of reference rotating with a uniform angular velocity f2j. The full SGS stress 
tensor (^) transforms in such a way that the SGS force, i.e. drij/dxj, in an inertial and in 
a rotating frame are the same. Horiuti fl^ has recently analyzed several SGS models and 



showed that some of them do not satify this condition. 

The velocity in a rotating frame, denoted by the asterisk, and the velocity in an inertial 
frame are related as follows 

Ui = -\- eimn^m^n • (^O) 

Transforming the substantial derivative to a rotating frame and using (|30|) gives 
d d d ^ d 

where the last two terms are the Coriolis force and the centrifugal force, respectively. 
Transforming the NS— a equation (^) to a rotating frame gives 



n^xi (31) 



d d dP d"^ d 

J "^333 



where r^^ is the expression (^) for 'star' quantities and the additional correction term 
appears 

eimn^mXlV'^V* - ekmj^m-^ - ekli^i-^ j . (33) 

The divergence of this correction stress tensor does not vanish; in fact. 

However, the first term in ( p^ combines with the Coriolis force in Eq. ( ^2] ) and the second 
term combines with the centrifugal force to yield the following form of the NS— a equation 
in the rotating frame: 

d d dU d 

J '' J J J 

where the total pressure H is 



LES equations with nonlinear dispersion Domaradzki & Holm February 6, 2008 9 



n = P 



a e 



dvl 
dx* 



(36) 



Remarks. This slight pressure shift is of no importance for incompressible fluid dynamics. 
Had we replaced ||Vvp in the second term of the kinetic energy norm (|TB|) by trace (e ■ e), 
as mentioned earlier, we would have removed this pressure shift at the cost of making the 
motion equation and stress tensor slightly more complicated. We shall decline this option 
here; however, it shall be required in the compressible case. 

Note that in the transformed motion equation (p5| ) the Coriolis force is computed using 
filtered velocity v*. This is the physical velocity of the fluid parcels in the rotating frame. 

One may also obtain the nondissipative part of the NS— a motion equation written in a 
rotating frame as in (^) by first substituting Vi = v* + eimn^mXn into Hamilton's principle 
in which the Lagrangian is given by the kinetic energy (|18D and then taking variations in 

0- 



the Euler-Poincare framework, as discussed in 



5 Relation to generalized similarity models 

We shall discuss the relation of the NS— a model to other SGS models in terms of the in- 



version of the filtering operation (0). In generalized similarity models for the NS equations 
the full NS velocity Ui in the expression for the SGS stress tensor (^) is obtained by an 
inversion of the filtering operation (0), sometimes called deconvolution. The deconvolu- 
tion is frequently encountered in image processing and its properties are discussed in many 
textbooks on this subject, e.g. Il^, 11^. In general, the inversion of the continuous filter is 



a ill-conditioned and often a singular problem precluding a unique solution. This difficulty 
leads to a concept of an approximate or regularized deconvolution such as a pseudo-inverse 
through a singular value decomposition approach . In the context of subgrid-scale model- 
ing the concept of deconvolution was discussed previously by several authors. Germano pO| 



introduced an exponential filter for which the filtering operation can be inverted exactly, 
i.e. the function f{x) in (|21| ) can be represented as a closed form expression in terms of 
the filtered function f{x). Using this representation one obtains a formal expression for the 
subgrid-scale stress tensor as a function of the filtered field. Germano shows that for 
this filter the dominant term in the expression for the subgrid scale stress has the form of a 
nonlinear velocity gradient model 



It] 



dui duj 
dxk dxk 



(37) 



This form, considered first by Leonard 
al. 



has been more recently investigated by Liu et 
and Winckelmans et al. |2^. In addition 
to the exact deconvolution procedure of Germano |2y] and Leonard other approaches 



22 1, Borue and Orszag |23], Leonard [24 
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have been proposed and utilized recently. Geurts ||2^ constructs an approximate inverse 
operator by requiring that polynomials up to a certain order are recovered exactly from 
their filtered counterparts. Stolz and Adams |27| use a formal power series expansion of the 
inverse operator and obtain the approximate inverse by truncating the expansion. Similar 
methods are used by Horiuti leading to a multi-level filtered model. The deconvolution 
is also an integral part of the SGS estimation model of Domaradzki et al. |3^. In the 

spectral space version of the estimation model an exact deconvolution is possible In 
the physical space representation the approximate deconvolution is accomplished by solving 
a sequence of tridiagonal equations 



Here, following Stolz and Adams and Horiuti IT^ we shall consider a formal inverse 



of the linear operator Lq in (|2T|) as a power series expansion 



La' = (/-(/- Lg))-' = I+{I~Lg) + {I- Lg) 

where / is the identity operator. To first order in the expansion 
approximated as 



+ ... (38) 
'^) the velocity Ui is 



Ui 



{ui - Ui) . (39) 

The second term in the approximation (^) is obtained by filtering of the difference {ui — Ui) 
which can be written using the Taylor expansion as 



Ui 



Ui 




X + x ) ) G(x' ) d^x' 



1 d'^Uj 

2 dxudxr 



+ . . . 



(40) 
(41) 



where the derivatives are evaluated at x and the above relation is valid for both the symmetric 
top hat filter or the Gaussian filter with the width A in each Cartesian direction. Using this 
relation the approximate inversion (|39|) becomes 



Ui 



Ui 



24^^^ 



(42) 



upon again invoking commutation of filtering and differentiation. By comparing (ff3) for the 



approximate inversion and (plf ) for the definition of Vi in the NS— a model, we recognize that 
Vi and Vi are related by the approximate inverse (^) with = A^/24. We shall pursue this 
relation farther by comparing the NS— a model with LES generalized similarity models. 

Employing the approximate inverse (^2|) in the expression for the SGS stress tensor ( ^2] ) 
one finds 
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The generahzed similarity tensor r/j'" in expression (^3|) is symmetric and Galilean invariant. 
On comparing the last equation with the expression for r/^" in Eq. (PP]), we find as the most 
prominent difference the product derivative term in the filtered alpha-model equation 



dxj dxi 

where r^p"^ is expression ( |43| ) written in terms of Vi and the 'hat' filter, instead of Ui and 
the 'bar' filter. Also f^^ is the 'hat '-filtered version of r^" in Eq. (26). According to Eq. 



(0) the SGS stress tensor for the filtered NS— a equation (^) is the sum of the similarity 
stress tensor for NS and the filtered SGS stress tensor (^) for the original NS— a equation 
(P5|). The last equality is the Germano identity |^ and follows from the the definitions 
ofr---, rf^, andr/;". 

We note that the important new term appearing in all SGS stress tensor expressions for 
the alpha-model, namely the product of derivatives, is different from the similar product of 
derivatives in the nonlinear models ([37|). 

The NS-a model is reminiscent, but not equivalent to the generalized similarity models. 
This is encouraging because it may avoid difficulties encountered by the pure similarity 
models. One of the most serious problems encountered by the pure similarity models is 
insufficient SGS dissipation caused by neglecting in the modeling scales which cannot be 
resolved on a given LES grid. This fact has been long recognized, e.g. by Zhou et al. [B2 



and more recently by Domaradzki and Loh . Indeed, consider the most favorable situation 



when the filtering operation can be inverted exactly, i.e. the function f{x) in (pl|) can be 
represented as a closed form expression in terms of the filtered function f{x). In such a 
case one immediately obtains an expression for the subgrid-scale stress tensor as a function 
of the filtered field from the definition (^). However, the apparent simplicity of such an 
exact approach to subgrid-scale modeling is deceptive. This is because in the procedure 
f{x) and f{x) have the same spectral support since the same LES mesh is used to represent 
both functions. Therefore, no modeling of the subgrid scales is actually performed in the 
deconvolution procedure. As a result the similarity SGS stress tensor is unable to account 
for the dissipative effects of the nonlinear interactions between the resolved scales and the 
unknown, subgrid scales which cannot be represented on a given LES mesh. This implies 
that a good SGS model should either explicitly model subgrid scales or possess extra terms 
in addition to the similarity part. The former approach is used in the SGS estimation 



modeling ||2^, ^] and the additional terms are encountered in the mixed models |33 



35 as well as in the NS-a model 



Section |5| summary. We have identified the 'hat '-filtered NS— a equations with a class 
of generalized LES similarity models in which the SGS stress tensor possesses additional 
terms involving products of derivatives, as in t;". Applying the approximate- inverse meth- 
ods of deconvolution to this class of models restores the Kelvin circulation theorem to the 
deconvolved solution. The corresponding dynamics of the deconvolved solution is the NS— a 



model, whose analytical and numerical properties are summarized in |T0 
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6 Spectral space interpretation 

The nonlinear term in the NS— a equation (^) transforms to spectral (Fourier) space as 

F.T. -^vi{x.)vj{x.) =ikj J vi{p)vj(k-p)d^p = ikj J vi(k-p)vj{p)d^p. (45) 

Here the spectral representation is indicated by the dependence on wavenumbers. The 
incompressibility condition requires projection of (^3|) onto a plane normal to the vector k 
which is accomplished by applying the tensor 



Pim{k) = {6im - kikm/k^). (46) 



In spectral space one has 



Vj (k) = (k) + a^kHj (k) , (47) 
so the projection of the nonlinear term can be written as 

ii /p^^.,„(k - p) + J i'p^ip)^^^§^'j . (48) 

This expression allows us to interpret the alpha-model in terms of triad interactions. 
Qualitatively, the factor 1/(1 + a^/c^) attenuates amplitudes of modes with wavenumbers 
k > 1/a and leaves modes with k <l/a largely unaffected. 

Consider first a mode with k < 1/a. Its triadic interactions with all other modes p and 
q = k — p for which p,q < 1/a are only slightly affected by the attenuation factors and so 
are essentially the same as for the Navier-Stokes dynamics. If one of the wavenumbers p, q 
is greater than 1/a its effect on the mode k is weakened by one of the attenuation factors 
in (^). In particular, the highly nonlocal interactions of the eddy viscosity type, k < 1/a 
and p,q ^ 1/a, are removed from the dynamics. 

Consider now a mode with > 1/a. In this case all local triads with p,q ^ 0{k) are 
attenuated and the nonlocal triads with p,q ^ 1/a are again removed from the dynamics. 
However, not all nonlocal triads affecting the mode k are neglected. Indeed, for k > 1/a it 
is possible to have triads with p < 1/a < k and q ~ 0{k) and then at least one of the terms 
in ( ^8]) is not directly affected by the attenuation factors. The NS— a model thus retains 
nonlocal sweeping effects of large scales (either p, or q < 1/a < k) acting on small scales 
k > 1/a. 



Section ^ summary. The NS— a model is equivalent to the Navier-Stokes dynamics for 
all triads with k,p,q < 1/a, it consistently neglects the nonlocal interactions of the eddy 
viscosity type for modes p,q > 1/a, and retains nonlocal sweeping effects of large eddies on 
small eddies {poTq<l/a<k). 
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